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ABSTRACT 

We present spectroscopic observations of the remaining four candidate white dwarfs in Praesepe. All 
four candidates are quasars with redshifts between 0.8 and 2.8. One quasar, LB 6072, is observed to 
have a strong metal- line absorption system blueward of the quasar redshift. The lack of additional 
white dwarfs in Praesepe leaves the total known white-dwarf population of the cluster at five, well be- 
low the number expected from commonly-assumed initial mass functions, though several undiscovered 
cluster WDs may lie in the outer regions of the cluster. All known Praesepe member white dwarfs are 
concentrated within 0?6 of the cluster center, and the radial profile of cluster white dwarfs is quite 
similar to the profile of massive cluster stars. This profile is mildly inconsistent with that of ~ IM© 
cluster stars and suggests that the white dwarfs did not receive a velocity kick during the progenitor 
star's mass loss phases. If complete, the observed Praesepe white dwarf population is consistent with 
a steeper high-end initial-mass function than commonly assumed, though the calculated slopes are in- 
consistent with the present-day mass function of Praesepe. Searches for white dwarfs outside the core 
of Praesepe and further study of the white dwarf populations of additional open clusters is necessary 
to constrain further the underlying cause of the white dwarf deficit. 

Subject headings: white dwarfs — open clusters and associations: individual (Praesepe) — quasars: 
general — quasars: absorption lines — stars: luminosity function, mass function 



1. INTRODUCTION 

White dwarfs (WDs) are the final stage of stellar evo- 
lution for the vast majority of stars. WDs are known 
to exist in open clusters with ma in-sequence turnoff 
masses of ~ 5Mq (e.g. NGC 2516. lKoester fc ReimersI 
119961) and may have progenitors with zero-age main se- 
quence masses as large as 8Mq. However, the WD 
population of the Hyades shows a deficit of WDs com- 
pared with the number that would be expected given 
the present-day mass function and reasonable assump- 
tions for the shape of t he initial-mass function (IMF) 
IjWeidemann et al] ll992"l. and no Hyades WDs have in- 
ferred progenitor masses > 3.6AfQ (Claver et al. 2001, 
hereafter CLBKOl). 

There arc at least three explanations for this white 
dwarf deficit. First, the "missing" WDs may be hidden 
in unresolved binary star systems. A calculation of the 
number of hidden WDs finds that this can explain many, 
but n ot all, of the missing WDs (Williams 2004, hereafter 
lW04fl . Second, dynamical evolution of the open clus- 
ter m ay remove WDs from the cluster (|Weidemann et aLl 
Il992t) . This effect is enhanced if the WDs receive a veloc- 
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ity " kick" during the mas s-loss phase of the progenitor 
star IjFellhauer et alJl200.'J) . Third, the higher-mass end 
of the initial-mass function (IMF) in the Hyades may 
have been steeper than commonly-observed IMFs, lack- 
ing stars with M > AMq. ^^_^^^ 

In part to address this issue, ICLBKOJ studied the 
WD population of the Praesepe open cluster. Praesepe 
has an age (^ 625 Myr) and metalhcity ( Z = 0.024 ) 
indistinguishable from that of the Hyades l)CLBK01() . 
Praesepe has a total mass of ^ 600M(7^ and a half- 
mass radius ~ 3.9 pc l|Adams et all 1200^ . more mas- 
sive and more compact than the Hyades, which has a 
tot al mass ^ 4:00Mp and a half-mass radius of « 5.7 
pc IjPerrvman et al.lll998() . so losses of WDs due to dy- 
namical evolution should be less in Praesepe than in the 
Hyades. 

In Praesepe, iCLBKOll found one WD with a massive 
progenitor, LB 1847 = EG 60, M^ = 4.17Af0, and four 
WDs with less-massive progenitors. Four other Prae- 
sepe WD candidates were not observed spectroscopically 
in that study due to their faintness; two of these (LB 
1839 and LB 6072) have photometric properties consis- 
tent with cluster membership, a large WD mass and a 
progenitor mass > 5Mq. If these objects are indeed clus- 
ter member WDs, they would be useful objects for ex- 
ploring the value of Merit and would mostly resolve the 
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Table 1. Observing log. 



WD 0834+209 



UT Date Telescope Object Exp time (s 



08 Dec 2002 Keck 



12 Feb 2004 
22 Mar 2004 



Keck 
Magellan 



WD 0834+209 
LB 1839 
LB 6072 
LB 6037 



600 
600 
600 
600 



Praesepe WD deficit at these higher masses l|W04() . 

As part of our ongoing study of open cluster WDs, we 
have obtained spectra of these four objects. §2 details 
the observations and spectroscopic identifications of each 
object, and §3 discusses the implications of this study. 

2. OBSERVATIONS AND DATA REDUCTION 

Spectra of Praesepe WD candidates were obtained be- 
tween 2003 and 2004 with the Keck I 10m telescope and 
the Magellan Baade 6.5m telescope. An observing log is 
presented in Tabled] We note that the coordinates for 
LB 6037, LB 6 072, LB 1839, and LB 1876 given in Table 
3 of lCLBKOH are labeled as J2000 coordinates, though 
the B1950 coordinates are published. For th e sake of 
clarity, J2000 coordinates for all the objects in lCLBKOll 
are given in Table [3 

The Keck observations u sed the upgrade d blue camera 
of the LRIS spectrograph IjOke et al-lllOOSJ) with the 400 
1 mm~^, 3400A blaze grism and a l"-wide longslit at 
parallactic angle. The Magellan observation used the 
IMACS spectrograph with the f/2 camera, the 200 1 
mm~^, 66OOA blaze grism, and the l"-wide slit in slit- 
viewing mode at parallactic angle. The resulting spec- 
tral resolution was 6A for the Keck data and 7A for the 
Magellan data. The spectra were reduced using standard 
NOAO IRAF routines. Relative flux calibration was ob- 
tained by applying the response function obtained from 
observations of spectrophotometric standards obtained 
on the corresponding observing run. 

Spectra for each of the four remaining white-dwarf can- 
didates are presented in Fig. ^ As can be seen, none of 
the four objects are WDs; all are QSOs. Redshifts for 
the QSOs were determined by cross-correlating the ob- 
ject spectrum with the SDSS composite QSO spectrum 
of Vanden Berk ct al. ( 2001) using the fxcor package in 
IRAF. For illustrative purposes, an example of the cross- 
correlation values for LB 6072 are shown Fig. [3 The 
spectral identifications and redshifts of each object are 
summarized in Table |5| Each of these quasars is also 
assigned a QSO catalog designation (listed in the first 
column of the table). 

Special difficulties were posed by LB 6 037. The co- 
ordi nates for this object in SIMBAD and iCLBKOll are 
from i Wagner et aO l|1986() : however, the object at these 
coordinates is not a blue object but a K star. In order 
to find any nearby blue objects, we blinked the blue and 
red second-generation Digitized Sky Survey plates from 
the ESO Onhne DSS server. A point source ~ 1' west 
of the given coordinates stands out as being far bluer 
than anything else in the field. A spectrum of the blue 
source obtained at Magellan reveals it to be a QSO. We 
therefore identify this blue object as LB 6037 and give 
these coordinates in Table |5] A finding chart is pre- 








Fig. 1. — Spectra of the four candidate Praesepe white dwarfs. 
Spectra are unbinnod and in (relative) units of f\ . 




Fig. 2. — Correlation values for LB 6072. The correlation value 
as a function of redshift for cross-corre lation of LB 6072 with the 
SDSS composite quasar spectrum from lVanden Berk et al.i (i2001ll . 
The sharp, high peak is the QSO redshift, 2 = 1.975. 



sented in Fig. O Due to the lower signal-to-noise of this 
spectrum, the cross-correlation routine failed to converge 
satisfactorily. Manual line identifications give a redshift 
of z = 2.78. 

LB 1839 was identified by ICLBKOll as a proper mo- 
tion member of Praesepe, though it is actually a QSO. 
As quasars have zero proper motion, this is an obvious 
error in the proper motion measurement. Measurements 
of "significant" proper motions for QSOs have been noted 
in a cross-correl ation of the USNO-B and SDSS catalogs 
(jGould fc Kollm eier 2003) ^'^d merely illustrate the diffi- 
culty in measuring small proper motions of faint objects. 

3. DISCUSSION 
3.1. Praesepe WD candidates and sample completeness 
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Table 2. Spectroscopic Identifications 



Object 



Other Names RA {J2000) Dec (J2000) V 



B-V^ 



ID 



QSO J0837+2040 


WD 0834+209'' 


8 37 11.5 


20 40 57.1 


18.83 


0.34 


QSO 


1.0916 + 0.0013 


QSO J0838+1933 


LB 1839 


8 38 37.2 


19 33 13.2 


18.83 


0.23 


QSO 


0.8622 ± 0.0009 


QSO J0843+1955 


LB 6037 


8 43 13.3 


19 55 59 


18.98 


0.37 


QSO 


2.78 + 0.01 


QSO J0843+1937 


LB 6072 


8 43 45.7 


19 37 24.2 


18.73 


0.24 


QSO 


1.9749 ±0.0004^= 


WD 0836+197 


LB 5893 


8 39 36.5 


19 30 43.2 


17.63 


-0.06 


WD 




WD 0836+201 


LB 393, EG 61 


8 39 45.6 


20 00 15.6 


17.96 


0.12 


WD 




WD 0836+199 


LB 1847, EG 60 


8 39 47.2 


19 46 11.8 


18.33 


0.08 


WD 




WD 0837+199 


LB 390, EG 59 


8 40 28.1 


19 43 34.4 


17.60 


0.04 


WD 




WD 0838+202 




8 41 39.8 


19 00 07.5 


17.96 


0.20 


WD 




WD 0840+200'* 


LB 1876 


8 42 52 


19 51 12 


17.69 


0.15 


WD 





''As given in CLBKOy, should not be used 

^Icr errors calculated in fxcor 

"^ Strong metal absorption Une system at 2 = 1.82 

•^Incorrectly Usted in fCLBKOll Table 3 as "WD 0837+202" 




Fig. 3. — Finder chart for LB 6037, indicated by vertical tick 
marks. The arrow indicates a K star at the coordinates given for 
LB 6037 in Wagner ct al. 1 1983). The image is a 2' X 2' field taken 
from the blue plate of the second-generation Digitized Sky Survey. 
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Fig. 4. — Mosaicked Digitized Sky Survey imag e of Prae sepe. 
Crosses indicate locations of Praesepe WDs from IGLBKOH the 
square indicates the CLBKOl non-Praesepe WD; open circles mark 
the location of the QSOs from this stu dy. The outlined region 
represents the approximate coverage of ICLBKOll The image is 
~ 2° on a side. 



Any conclusions about the Praesepe WD population 
will depend on the completeness of the WD candidate 
lists. CLBKOl do not discuss the completeness of their 
sample. They also included a few previously known WD 
candidates outside their survey area (outlined in Figure 
4). However, the CLBKOl photometry is complete to 
V ^ 19.5, more than 1 mag fainter than the faintest 
Praesepe WD, so it is not unreasonable to assume that 
the sample is complete (aside from the complication of 
WDs in binary systems) in the central 2.1 square degrees 
(^ 50' radius) of Praesepe. 

3.2. The Deficit of Praesepe WDs 

Memb ers of Praesepe are k nown to exist out to a radius 
of ~ 8° IjAdams et al.ll200'2l) . so it is quite possible that 
more cluster WDs exist in the outer regions of the cluster, 
just as the WD population in the open clust er M67 is ob- 
serve d to extend to the largest cluster radii (Richer et al.l 
I1998J) . The potential number of Praesepe WDs outside 
the surveyed area can be estimated by several means. Us- 



ing the p arameters of the bes t-fitting King model deter- 
mined bv lAdams et all lj2002(l . we find that ~ 85% of the 
stellar content of Praesepe is outs ide the central 50'. Us- 
ing the profile of "^ IMq stars in lRaboud fc MermilliodI 
(|1998b^, we estimate that ~ 80% of the stellar content is 
outside the surveyed regions. Therefore, if the WD pop- 
ulation of Praesepe has the same radial distribution as 
that of the majority of the stellar population (discussed 
in ^'S.'Iti . Praesepe may have a total population of ^ 25 
to 30 WDs. 

The number of expec ted, observable WDs in Praesepe 
was calculated in lW04l to be seven to 2 1 WDs. This cal- 
culatio n used the luminosity function of lJones fc Stauffen 
ljl991(l . which covered the central ~ 16 sq. degrees 
of Praesepe. The large range in the number of ex- 
pected WDs arises primarily from the assumed IMF. 
If this calculation is scaled by the radial profile of 
iRaboud &: MermilliodI l)1998bf) under the assumption 
that the WD radial profile follows that of the ~ IMq 
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stars, we estimate that ^^ 1 to 5 WDs should be found 
in the central region studied by;CLBKOH, in agreement 
with the five known Praesepe WDs. If, however, these 
five known Praesepe WDs represent the entire WD popu- 
lation of the cluster, then Praesepe suffers from the same 
deficit of WDs as the Hyades (Weid emann et al.,!:1992il . 

WD Velocity Kicks 

One explanation for the deficit of WDs in open clus- 
ters is that during the mass-loss phase of the pro- 
genitor star's stellar evolution, the future WD is ac- 
celerated due to asymmet ric mass loss, resulting in a 
small net velocity "kick" IjWeidemann et al1ll992t) . A 
net velocity of just a few km/s is sufficient to remove 
the majority of WDs from an op en cluster on short 
time scales (jFellhauer et al.l [200,'^ . In N-body simu- 
lations not considering velocity kicks, the WD popu- 
lations of open clusters with ages similar to Praesepe 
have radial profiles very si milar to those of the lumi- 
nous stars and subgiants (' Portegies Zwart et al.1 l200lt 
iBaumeardtTM akino 2003). 

Figure 0| shows the location of the five Praesepe WD 
candidates with respect to the cluster. The figure shows 
a 2° X 2° region centered on Praese pe, which agrees we ll 
with the cluster core radius of ~ 1° (jAdams et al.ll200^ . 
Also indicated on the figure is the region of sky studied 
by CLBKO^ All of the cluster member WDs are located 
within ^ 35' of the cluster center, even though the region 
studied in CLBKOl extends > 50' out from the center. 

Figure [Sji compares the cumulative distribution of 
WDs as a function of radius from the cluster center to 
th e distribution of cluster stars in terior to 50' from data 
bv lRaboud fc MermilliodI l|1998bf) . From the figure it ap- 
pears that the WD distribution follows that of the mas- 
sive cluster stars, exactly as predicted for a WD popu- 
lation with no net velocity kick. A K-S test comparing 
the WD radial distribution to that of the stars finds that 
the WD distribution is consistent with that of the stars 
with masses > 1.6Mq and weakly inconsistent with the 
distribution of lower mass stars. 

To test further the significance of the WD radial 
distribution, two Monte Carlo simulation were per- 
formed. In the first, five stars were drawn from the 
central 50' radial distr i bution of the < 1.2Af0 stars in 
IRaboud fc MermilliodI l|1998hD . Out of 10,000 realiza- 
tions, 13% had all five selected stars interior to the 
35' extent of the known WD population. The sec- 
ond Monte Carlo simulation drew stars from the ra- 
dial distribution of M < 1.2 Mf:^ stars from the entire 
IRaboud fc MermilliodI flOOSb") sample until five stars in- 
terior to 35' were found. Of these, only 15% of the re- 
alizations were found to have no additional stars in the 
annulus between 35' and 50'. 

These simulations hint that the WD population is more 
consistent with the giants and subgiants of Praesepe than 
with the lower mass, as-yet unevolved cluster members, 
as predicted by the N-body calculations that do not in- 
clude a velocity kick. Figure|SlD shows the resulting radial 
distribution of WDs compared with those for other Prae- 
sepe stars if the five known WDs represent the entire WD 
population of Praesepe. In this case, the similarity of the 
WDs with subgiant and giant Praesepe members is even 
more striking. 

The concentration of WDs toward the center of the 



cluster may not be inconsistent with the concept of ve- 
locity kicks if all the WDs with non-zero kicks have al- 
ready been lost from the cluster (or are at larger radii). 
It would be useful for N-body simulations of open clus- 
ters to explore how the apparent radial distribution of 
WDs varies with the strength of the applied kick. 

We also note that the Pleiad WD (LB 1497) is lo- 
cated ~ 1?4 away from the tightly-concentrated core 
of the cluster, well outside the half-mass radius (w 
53') but w ithin the Pleiades tidal radius (« 7?4) 
(jRaboud fc Mermilliod 1998a) . This contrasts with the 
central concentration of WDs observed in Praesepe and 
further serves to muddy the waters. 

Shape of the Initial-Mass Function 

It is also possible that the dearth of WDs in Praesepe 
could be due to a steeper IMF than normally consi dered. 
To explore this possibility, we use the code from IW04I 
to calculate the WD population of Praesepe for various 
IMF slopes. The parameters used to normalize the IMF 
and the WD detection criter ia are the same as those for 
Praesepe presented in |W0_-|: Merit is taken to be SMq, 
and the binary fraction is taken to be 0.4. Two differ- 
ent IMFs were explored - a power-law IMF of the form 
$(M) ex M~ ^ ^+'"^ a nd the broken power-law IMF from 
INavlor et al.l lj2002() . with a steep high-mass power law 
IMF (F = Fi for M > IMq) and a flat low-mass slope 
(F = F2 = 0.2 for M < IMq). 

The Praesepe member W D 0837-1-199 has a progenitor 
mass < 4.17Mq ||CLBK01D . We therefore calculate the 
likelihood that, for a given IMF slope F, no WDs with 
progenitor masses > 4.17M0 are observed. The value of 
F was then varied to find the slopes corresponding to con- 
fidence levels of 60%, 90%, and 95%. Results are given 
in Table |21 The calculations show that the lack of high- 
progenitor-mass WDs, the lack of WDs observed in bi- 
nary systems, and the overall deficit of WDs in Praesepe 
can be explained by a steep high-mass IMF. Assuming 
Poisson statistics, both the F = 2.55 power-law IMF and 
the Fi — 2.75 broken power-law IMF are remarkably con- 
sistent with the Praesepe WD population, though such 
steep IMFs ar e contrary to what is observed in young er 
open clusters llKroupall200l IKroupa fc Weidne r 2003). 

The present-day mass function (PDMF ) of Praesepe 
has b een discussed by several authors. lAdams et alJ 
l|200l find a PDMF slope of F = 0.6 for IMq > M > 
O.4M0, somewhat shallower than the Salpeter slope of 
F = 1.35 and far shallower than the IMFs needed to ex- 
plain the WD deficit. Williams. Ricke. fc Stau ffcr ( 1995|) 
find a steeper IMF for stars more massive than 0.6A/q, 
with F = 1.7, though this is still shallower than the IMF 
required by our calculations. 

3.3. Future work 

In order to determine the significance of the observed 
WD deficit in Praesepe, and therefore determine the va- 
lidity of the above discussion, it is necessary to conduct 
a search for Praesepe WDs out to larger radii. Including 
proper motions in such a study would be very helpful 
in reducing the large number of background and extra- 
galactic objects. Unfortunately, existing proper motion 
surv eys of Praesepe do no t extend quite faint enough 
(e.g. I.Tones fc Staufe^[T99l . 
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Fig. 5. — Radial distribution of stars and WDs in Praesepe (a) interior to 50' and (b) over the entire cluster. Large, filled circles show 
the cumulative radial profile of the five known WDs. The remaining data are taken from Raboud & Mermilliod 1 1998b), show the profile 
of stars with M > 2.5Mq (triangles), 2.5Mq > M > 1.6Mq (crosses), l.GM© > M > IMq (open squares), M < IMq (filled squares). 
The zero-age main-sequence masses of Praesepe giants are ^ 2.7Mq. 



Table 3. Calculation of the Praesepe WD population 



NwD 













Observed 




Ibtal 


O bserved 


IMF'' 


Confidence 


r 


Total 


Observed 


in Binaries 


{Ah 


> 4.17M0) 


(Mi > 4.17 Mq) 


P 


60% 


2.55 


5.7 


4.1 


1.6 




1.9 


0.95 


P 


90% 


2.02 


12.1 


8.0 


2.8 




4.8 


2.3 


P 


95% 


1.87 


15.0 


9.7 


3.3 




6.2 


3.0 


N 


60% 


2.75 


6.4 


3.4 


0.7 




2.1 


0.9 


N 


90% 


2.16 


13.6 


7.2 


1.3 




5.1 


2.3 


N 


95% 


2.75 


16.7 


8.8 


1.6 




6.6 


3.0 



ip= Power-Law IMF, N = INavlor et"an <2002) IMF 



Examination of the complete WD content of additional 
clusters should be able to constrain these hypotheses 
even further. Comparison of the WD populations of 
similar-age clusters of varying total masses can deter- 
mine the degree of any dynamical contributions to the 
WD deficit, as more massive clusters should retain a 
larger fraction of WDs if dynamical effects dominate the 
WD deficit. Examination of young stellar clusters with 
masses similar to Praesepe and the Hyades may help to 
constrain the high-mass IMFs, especially if a steep IMF 
is observed in these clusters. 
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